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Abstract 

A series of aerodynamic heating tests was conducted 
on a 70-deg sphere-cone planetary entry vehicle 
model in a Mach 10 perfect-gas wind tunnel at 
freestream Reynolds numbers based on diameter of 
8.23xl0 4 to 3.15x10 s . Surface heating distributions 
were determined from temperature time-histories 
measured on the model and on its support sting using 
thin-film resistance gages. The experimental heating 
data were compared to computations made using an 
axisymmetric/2D, laminar, perfect-gas Navier-Stokes 
solver. Agreement between computational and 
experimental heating distributions to within, or slightly 
greater than, the experimental uncertainty was obtained 
on the forebody and afterbody of the entry vehicle as 
well as on the sting upstream of the free-shear-layer 
reattachment point. However, the distributions began 
to diverge near the reattachment point, with the 
experimental heating becoming increasingly greater than 
the computed heating with distance downstream from 
the reattachment point. It was concluded that this 
divergence was due to transition of the wake free shear 
layer just upstream of the reattachment point on the 
sting. 

Nomenclature 

B c H heating bias error 

C H Stanton number, qt [p^fhg—h)] 

d molecular diameter (m) 

h enthalpy (J/kg) 

Kn Knudsen number 

k thermal conductivity (W/m-K) 

M Mach number 
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N number density (1/m 3 ) 

P c H heating precision error 

p pressure (N/m 2 ) 

Q heat (J/m 2 ) 

q heat transfer rate (W/m 2 ) 

R radius (m) 

Re Reynolds number 

S distance along model surface (m) 

T temperature (K) 

t time (sec) 

U ^ freestream velocity (m/sec) 

a thermal diffusivity (m 2 /sec) 

P thermal product, al 4k: (W-sec 1/2 /m 2 -K) 

AC h heating uncertainty 

X correction factor (1/K) 

X mfp mean free path (m) 

p density (kg/m 3 ) 

Subscripts: 

0 stagnation 

1 freestream 

2 post-shock 

D diameter 

w wall 

Introduction 

Wake flow behavior is one of the important 
issues which must be considered in the design of 
planetary entry vehicles such as the Mars Pathfinder 
probe 1 . The nature of the wake flow dictates payload 
size, placement and shielding requirements. Accurate 
characterization of the wake flow behavior is becoming 
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more important due to the increasing reliance in 
planetary mission planning on aerobraking, which 
produces much more severe aerothermodynamic loads on 
an entry vehicle. 

In this study, the effects of transition in the 
wake of a blunt body were investigated through heat- 
transfer testing of a 70-deg sphere-cone entry vehicle 
configuration in a perfect-gas hypersonic wind tunnel. 
Other recent studies 2 ' 15 have also addressed the effects of 
nonequilibrium thermochemistry, rarefaction, and 
transition to turbulence on wake flow behavior. 

Experiment Description 

Test Models 

The test model geometry was a 70-deg sphere- 
cone forebody of 1.0-in. radius with a 40-deg cone- 
frustrum afterbody (Fig. 1). The forebody corner radius 
was 0.05 forebody radii, and the radius of the 
hemispherical nose was 0.5 forebody radii. The 
afterbody frustrum radius was 0.6 forebody radii. 

A heat-transfer test model was fabricated from 
Macor, a thermally-insulative, machinable glass-ceramic 
material. The model was instrumented with fast- 
response (< 1 psec) thin-film temperature resistance 
gages. There were a total of 37 gages on the forebody, 
afterbody and base of the model (Fig. 2). The support 
sting for this model was fabricated from stainless steel. 
A lengthwise slot into which a contoured Macor insert 
with 33 additional thin-film gages was fitted was 
machined into the sting (Fig. 2). The support sting 
was fitted to a sting adapter and the facility support 
barrel at a distance of 5.25 forebody radii downstream of 
the model. An additional uninstrumented model and 
sting were fabricated from aluminum for use in oil-flow 
tests. 

Facility Description 

Aerothermodynamic testing was conducted in 
the NASA Langley Research Center (LaRC) 31-Inch 
Mach 10 Air Tunnel 16 (Fig. 3). This facility is a 
conventional perfect-gas hypersonic blowdown wind 
tunnel in which air is used as the test gas. The 31-Inch 
Mach 10 Air Tunnel has been calibrated for operation at 
reservoir pressures of 0.86 MPa to 10 MPa at a 
reservoir temperature of 1000 K. These reservoir 
conditions produce freestream Reynolds numbers of 
from 1.62xl0 6 nr 1 to 6.20xl0 6 m _1 at a nominal 
freestream Mach number of 10. Facility run times for 
the current series of heat-transfer tests were limited to 3- 
5 seconds in order to avoid violation of the semi-infinite 


heat transfer assumption 17 , although much longer run 
times are possible for aerodynamic testing. Freestream 
conditions for this study are listed in Table 1. The run- 
to-run repeatability of the facility flow conditions was 
investigated, and it was found that the freestream 
conditions were repeatable to within less than ±2%. 


Data Reduction 

Heat-transfer rate time-histories were computed 
from thin-film gage temperature time-histories using 
the 1DHEAT data reduction code 18 . Heating time- 
histories were determined from the temperature time- 
histories by two methods: the Kendall-Dixon method 19 , 
which is a closed-form analytical scheme; and a 
numerical solution of the one-dimensional heat- 
conduction equation. Both methods are based on the 
assumption of one-dimensional heat conduction within 
the model substrate 17 . 

In the Kendall-Dixon method, the total heat 
energy added as a function of time is first computed by: 


Q(0 


P j, (T, ~ T t ) + (T j _ l - T { ) ^ 

^ i = 1 V K ~ h + V f n ~ h-\ 


( 1 ) 


After the total heat added has been computed, the heat- 
transfer rate is computed from 20 : 


^ , d Q(tJ -22„-8 _ Q, ,-4 + e„ +4 + 2 Q„ + S 
C 1 ( 0 = — = 7777 “ ( 2 ) 


dt 


40 At 


It is assumed in the derivation of the Kendall- 
Dixon method that the thermal properties of the model 
remain constant. However, the thermal properties of 
Macor are functions of temperature. Curve fits to 
experimental thermal properties data 18 are given by: 

p= 2543 (kg/m 3 ) (3) 

k = 0.33889 + 7. 4682x10“ 3 T - 1.6118xl0“ 5 T 2 
+ 1.2376xl0' 8 T 3 (W/m-K) 

(4) 

a = 1.3003xl0“ 6 - 2.2523xl0“ 9 T + 

(5) 

1.8571xl0“ 12 T 2 (m 2 / sec) 


The results from the Kendall-Dixon method are 
adjusted to account for the temperature dependence 
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through an empirical correction given by 18 : 


q = q ,[l + A(r - 300)1 (6) 

L var L const \ w J x / 


where the correction factor X is given by: 

A = 7.380xl(T 4 - 4.604xl(T 7 (T w - 300 ) (7) 


This empirical correction can be avoided 
through the use of the numerical method, in which the 
variation of thermal properties with temperature is 
included in the formulation of the method. In this 
numerical method, the one-dimensional heat conduction 
equation: 


P°n 


d t_ 
dt 


— fit—] 

dx V dx ) 


(8) 


is solved through a time-implicit, finite-volume 
discretization given by: 


[, )] 


( r r n T n \ 

u \ D+i - ) 

- 1 ' 2 { x r x i-i) 


i + H2 (*,+ "*,) 


/ y w j 1 n — 1 


The Kendall-Dixon method was used for rapid 
initial analysis of the experimental data. The validity of 
the empirical correction in this method was later verified 
by re-reducing the data using the finite-volume method. 

Experimental Results 

The 70-deg sphere-cone heat- transfer model was 
tested in the 31 -Inch Mach 10 Air Tunnel at diameter- 
based freestream Reynolds numbers of 8.23xl0 4 , 
1.62xl0 5 ' and 3.15xl0 5 . The experimental data are 
reported in terms of the Stanton number which is 
defined by: 


C 


H 


q 

Pi u i[K - K) 


( 10 ) 


The Stanton number remains nominally constant during 
a test, whereas the heat-transfer rate decreases as the wall 
surface temperature and enthalpy increases. 

Stanton number distributions normalized by 


the measured stagnation point values from tests at each 
of the three Reynolds number test points are plotted in 
Figs. 4 and 5. The complete distributions are shown on 
a log scale in Fig. 4, while the details of the forebody 
and wake distributions are shown separately on linear 
scale plots in Fig. 5. The Reynolds number had no 
effect on the normalized forebody heating distribution, 
but had a strong influence on the heating distribution in 
the wake, where the normalized heating was an 
increasing function of Reynolds number. The peak 
heating on the model sting varied from 8% of the 
forebody stagnation point heating at the lowest 
Reynolds number to 15% at the highest Reynolds 
number. Relative peak sting heating rates of this level 
are consistent with turbulent wakes 6 - 714 , whereas 
relative laminar peak sting heating levels are in the 4- 
5% range 12-15 . 

Several surface oil flow tests were also carried 
out using the uninstrumented aluminum model. From 
inspection of oil flow photographs from these tests, it 
was possible to identify the boundary-layer separation 
point on the afterbody and the free-shear-layer 
reattachment point on the sting. A sample surface oil 
flow photograph is shown in Fig. 6. It was found that 
the reattachment point on the sting moved a small 
distance upstream toward the model as the Reynolds 
number was increased. This behavior is typical of 
increasing turbulence in the free shear layer. As the 
turbulence increases, momentum dissipation causes its 
thickness to increase, and the thicker turbulent free shear 
layer then comes into contact with the sting further 
upstream than a laminar free shear layer would. 


Experimental Uncertainty 

Two primary sources of experimental 
uncertainty were considered: precision error, P CH , due to 

run-to-run repeatability of flow conditions and test data, 
and bias error, B ( H , due to uncertainty in the thermal 
properties of the Macor models. From these two 
sources, the overall uncertainty, AC H , was estimated as: 


AC, 

C„ 


= t„ 


/ B 'F 

U C,H 


C 


( P C,H V 


H ) 


c 


H J 


(ID 


As per the AIAA Standards Committee 
recommendations 21 , the value of the 95th percentile 
point of the two-tailed Student’s “t” distribution was set 
to t 95 = 2, which is the limit for an infinite number of 
data points. Due to the high flow quality of the 31 -Inch 
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Mach 10 Air Tunnel, the precision error in run-to-run 
repeatability was only on the order of +1 - 2 %. From 
thermal properties calibration test data reported in Ref. 
18, the thermal properties bias error was estimated to be 
on the order of ±3%. The overall total uncertainty for 
heating in this study was thus estimated to be on the 
order of ±7 - 8%. 

Computations 

Code Description 

Flow field solutions for the 70-deg sphere -cone 
geometry were generated using the laminar, 
axisymmetric/2D Navier-Stokes solver NEQ2D 22 . A 
perfect gas thermochemical model was used to reproduce 
conditions equivalent to those of the 31 -Inch Mach 10 
Air Tunnel tests, although NEQ2D does also 
incorporate nonequilibrium thermochemical models. 
The governing equations are solved in NEQ2D through 
the line Gauss-Seidel method. Inviscid fluxes are 
represented by a modified form 23 of the Beam-Warming 
flux-splitting technique which minimizes the numerical 
dissipation produced in boundary layers by the original 
form. Viscous fluxes are represented by central second- 
order differences. 

Grid Generation, Adaptation and Resolution 

An initial grid of 125 streamwise by 90 
normal points for the sphere-cone geometry was created 
using an elliptical grid generation algorithm 24 . 
Orthogonality and initial cell spacing at the wall were 
enforced through Steger-Sorenson 25 boundary 
conditions. During computation of a solution, the grid 
was aligned with the bow shock and cells were clustered 
at the wall using an algebraic adaptation algorithm 
developed by Gnoffo 26 . 

It was found that the wake computations were 
extremely sensitive to the structure of the wake grid. 
This sensitivity is due to the fact the wake free shear 
layer and recirculation vortex are viscous-dominated 
regions. Grid resolution in the wake vortex and free 
shear layer must be of the same order as found in 
attached wall boundary layers in order to properly 
resolve the flow field gradients in these regions. 

In order to optimize the wake grid structure, 
the Volume Grid Manipulator 27 (VGM) code was used 
to cluster grid points within the recirculation vortex and 
to align the streamwise grid lines with the free shear 
layer (Fig. 7) after an initial flow field solution was 
computed. This adaptation produced significant changes 
in the structure of the wake flow field. As seen in Fig. 


8, the free-shear-layer reattachment point on the wake- 
adapted grid moved downstream from the computed 
location on the unadapted grid and so the recirculation 
vortex grew larger. This growth in the wake vortex size 
is due the fact that numerical dissipation produced by 
the lack of resolution within the wake vortex is 
minimized by the adaptation of the wake grid through 
VGM. 

In addition to the effects of the wake grid 
structure on the computed solutions, the effects of grid 
point density were also investigated. Flow field 
solutions were found to be relatively insensitive to 
streamwise grid resolution as compared to the 
sensitivity to normal grid resolution. Therefore, to 
investigate the effects of normal grid point resolution 
on the solution, grids with 125 streamwise points and 
45, 179, and 357 normal points were constructed from 
the baseline (125 x 90) point shear-layer-adapted grid. 

The forebody computations (Fig. 9a) were 
relatively insensitive to normal grid point resolution. 
The greatest sensitivity was observed around the nose 
and corner heating peaks, where the heating rates 
decreased by less than 5% between the 45 and 90 normal 
points grids, and by ~1% between the 90 and 179 point 
grids. 

Wake computations (Fig. 9b) were much more 
sensitive to grid resolution than forebody computations, 
as would expected given the much larger viscous region 
to be resolved in the wake. The sensitivity was greatest 
on the part of the sting bounded by the recirculation 
vortex, while the results on the portions of the 
afterbody and sting over which the flow was attached 
showed less sensitivity. On the portion of the sting 
within the wake vortex, heating rates decreased by more 
than 10% between the 45 and 90 normal point grids, by 
5 - 10% between the 90 and 179 point grids, and by less 
than 5% between the 179 and 357 point grids. 

From these grid resolution results it was 
concluded that 90 normal points were more than 
sufficient to model the forebody flow and the attached 
flow on the afterbody and on the sting downstream of 
the vortex. However, acceptable results were achieved 
within the recirculation vortex only on the highest 
resolution, 357 normal point grid. 

Computational Results 

Solutions were computed for each of the three 
Reynolds number operating points on shear-layer- 
adapted, (125 x 357) point grids using the NEQ2D code. 
Freestream flow conditions were taken from Table 1, 
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and a uniform 300 K wall temperature was specified. 
These computations revealed complex recirculating flow 
patterns in the near wake of the body. As shown in 
Fig. 10, a large main vortex was generated between the 
free shear layer and body. Smaller counter-rotating 
vortices were also produced at the base of the vehicle 
and on the afterbody. The size of the main vortex 
increased with Reynolds number, which is the 
theoretically-expected behavior in a laminar wake: as the 
Reynolds number is increased, more momentum is 
added to the shear layer, and thus the reattachment point 
is pushed further downstream from the separation point 
on the body. The size of the smaller vortices also 
depended on Reynolds number, and at the lowest 
Reynolds number the smaller afterbody vortex 
disappeared. 

Heat-transfer distributions normalized by the 
computed stagnation point values for each of the three 
Reynolds number test cases are plotted in Figs. 1 1 and 
12. The computed distributions are shown on a log 
scale plot in Fig. 11, while the details of the forebody 
and wake distributions are shown separately on linear 
scale plots in Figs. 12a and 12b. As in the 
experiments, the Reynolds number had no effect on the 
normalized forebody heating distribution (the slight 
spread in the forebody distributions is attributed to the 
use of the same boundary layer grid spacing for each of 
the three different Reynolds number cases), while in the 
wake, the normalized peak sting heating increased with 
Reynolds number. 

Comparison of Experiment and 
Computations 

Experimental and computational heating 
distributions for each of the three Reynolds numbers are 
compared in Figs. 13-15. The dimensional values for 
the experimental heating are extrapolated to the uniform 
300 K wall temperature of the computations from the 
Stanton number defined in Eq. (10). Uncertainty 
bounds shown in these figures are based on the 
uncertainty analysis discussed previously. 

Agreement between the two sets of 
distributions was very good on the forebody for all test 
cases. Computed heating rates were consistently higher 
than the experimental data, but remained within or just 
slightly outside of the estimated experimental 
uncertainty bounds. 

In contrast to the forebody results, 
computational and experimental wake heating 
distributions differed significantly. Fair agreement was 


achieved on the afterbody and on the sting close to the 
model; however, the two distributions begin to diverge 
just upstream of the point where the free shear layer 
reattaches to the sting. Furthermore, while computed 
peak sting heating rates varied with Reynolds number to 
from 6% to 8% of the stagnation point heating, the 
experimental values ranged from 8% to 15% of the 
stagnation point heating. It was theorized that the 
divergence of the experimental and computational 
results was due to the growth of turbulence in the 
reattached boundary layer on the sting, and that 
transition occurred just upstream of the free shear layer 
reattachment point (Fig. 16), which is where the 
heating distributions begin to differ. 

While transition had a large effect on the flow 
downstream of reattachment, it had little effect on the 
flow within the wake vortex since the favorable pressure 
gradient produced as the reversed flow within the vortex 
accelerates from rest at the reattachment stagnation 
point tends to dissipate turbulence ingested into the 
wake vortex. For this reason, the computed laminar 
heating distributions upstream of reattachment and the 
computed vortex sizes agree fairly well with the 
experimental data. 

The theory that the reattachment region is 
where transition occurs is based on studies by Lees 28 , 
Demetriades 29 , and Zeiberg 30 . In these studies, it was 
observed that at low Reynolds numbers, transition first 
occurs far downstream from the reattachment point (or 
from the neck of the wake in the case of an actual 
vehicle, which would not have a sting and reattachment 
point). As the Reynolds number is increased, the 
transition point moves upstream until it reaches the 
reattachment point (or neck). Because of the increasing 
pressure gradient which is caused as the flow is 
decelerated and turned at reattachment, the reattachment 
region acts as a trigger for transition, and so the 
transition point remains fixed here across a wide range 
of Reynolds numbers. Eventually, the Reynolds 
number is high enough for transition to occur in the 
free shear layer, and the transition point begins to move 
upstream within the free shear layer. 

Comparisons between the computed 
streamlines and the oil flow photographs support the 
hypothesis that transition occurs near the reattachment 
point. Although the laminar computations predict a 
small downstream movement of the separation point 
with Reynolds number and the oil flow photographs 
reveal a small upstream movement, the actual locations 
were quite close (Fig. 17), which suggests that 
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transition cannot be occurring very far upstream of 
reattachment. Furthermore, the computed and observed 
locations of the boundary layer separation on the 
afterbody of the model are also close, which would 
suggest that the attached boundary layer on the model is 
laminar, as does the good agreement between measured 
and predicted heating distributions on the forebody 
experiment. 

The transition of the free shear layer to 
turbulent flow and the location of the transition point 
have been presented as hypotheses. Transition appears 
to be the best explanation for the discrepancies between 
predicted and measured heating levels on the sting, 
which are much greater than the estimated experimental 
uncertainty. However, without flow field diagnostics 
such as hot-wire surveys, or schlieren or electron beam 
photography, which were not available in this study, 
definitive conclusions about the state of the free shear 
layer cannot be made. Nevertheless, two other possible 
causes of disagreement between experiment and 
computation can be eliminated. These are the 
computational grid structure and rarefaction effects in 
the wake. 

While it has been shown that the wake 
computations can be very sensitive to the grid structure, 
all grid adaptations and grid resolution increases 
performed in this work had the effect of lowering the 
heating on the sting and moving the reattachment point 
further downstream. That is, optimization of the wake 
grid further increased the differences between 
computational and experimental results. Therefore, the 
grid structure cannot be the cause of the differences 
between experiment and computation. 

However, in regard to grid structure, it is 
important to note that the magnitudes of the computed 
laminar heating rates on an original, unadapted grid (Fig 
9b.) are comparable to those of the turbulent data (Fig. 
13b). Since laminar computations on non-optimized 
grids may, incorrectly, be in qualitative agreement with 
turbulent data, comparisons with these computations 
could lead to incorrect conclusions about the state of the 
wake. An unpublished preliminary comparison with 
computational results from a non-optimized grid was 
one of the reasons that the experimental data discussed 
in the present paper were originally thought to result 
from a laminar wake 4 . While other issues such as flow 
quality, rarefaction effects, and thermochemical models 
in the CFD codes must also be investigated, the grid 
structure should also be considered as one of the 
possible reasons for the disagreement seen in wake 


heating comparisons such as Refs. 10 and 13. 

In order to confirm that the Navier-Stokes 
equations for continuum flow were applicable for the 
low density flow in the wake of the model, local 
Knudsen number contours were computed. The local 
Knudsen number is given by: 


„ ''■MFP I , 

K n = — H V Pl 


where the mean free path is: 


A 


MFP 


1 

spl nd 2 N 


( 12 ) 


(13) 


According to Bird 31 , the continuum 
assumption begins to break down for local Knudsen 
numbers above 0.1, and for values above 0.2 the 
Navier-Stokes equations are no longer valid. 

Local Knudsen number contours for the three 
cases are shown in Figs. 18-20. At the high and 
intermediate Reynolds number cases, the local Knudsen 
numbers reach 0.1 only in a thin region near the 
separation point, while the local Knudsen number 
reaches the 0.2 level only in a small area around the 
separation point for the lowest Reynolds number case. 
Therefore, the Navier-Stokes equations are considered to 
be appropriate for the conditions of this study 


Summary and Conclusions 

Experiments and computations have been 
performed for a 70-deg sphere-cone entry vehicle 
geometry at Mach 10. The experiments were conducted 
in a perfect-gas wind tunnel, and the computations were 
made using a laminar, axisymmetric/2D Navier-Stokes 
solver with a perfect-gas thermochemical model. 
Comparisons were made between heating distributions 
from the experiments and computations. Agreement of 
the heating distributions to within the experimental 
uncertainty was achieved on the forebody of the model. 
In the wake, the experimental heating data began to 
diverge from the computational results just upstream of 
the free shear layer reattachment point on the sting. 
Computed peak heating rates on the sting varied with 
Reynolds number from 6% to 8% of the forebody 
stagnation point heating, while the peak sting heating 
from the experiment varied from 8% to 15% of the 
stagnation point heating. From these comparisons, it 
was concluded that in the experiments transition to 
turbulence in the wake occurred in the vicinity of the 
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free-shear-layer reattachment point on the model support 
sting. 
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Table 1: 31-Inch Mach 10 Air Tunnel Test Conditions 



Re 1D = 8.23xl0 4 

Re 1D = 1.62x10 s 

Re 1D = 3.15xl0 5 

Pr (Pa) 

69.0 

130.6 

242.0 

T, (K) 

53.31 

52.45 

51.48 

Pi (kg/m 3 ) 

0.00451 

0.00868 

0.01646 

U , (m/sec) 

1416 

1422 

1425 

Mi 

9.675 

9.795 

9.928 

h 0 " IU 98 K (^J(k§) 

0.756 

0.764 

0.767 

P 0,2 

8383 

16280 

31000 

H 

© 

1008 

1015 

1017 
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Figure 1 . Entry Vehicle Model Geometry 




Figure 3. 31-Inch Mach 10 Air Tunnel 
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Figure 4. Reynolds Number Effects on Experimental 
Heating Distributions 
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Figure 5a. Reynolds Number Effects on Experimental 
Heating Distributions - Forebody Details 



Figure 5b. Reynolds Number Effects on Experimental 
Heating Distributions - Wake Details 



Figure 6. Sample Surface Oil Flow Patterns, 
R e i, D = 8.23xl0 4 
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Figure 7: Shock-Aligned Grids Before and After Wake Adaptation Using VGM, Repo = 8.23x 10 4 
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Figure 8. Effects of Grid Adaptation on Wake Flow, Rej D = 8.23xl0 4 
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Figure 9a. Grid Resolution Effects on Computed 
Forebody Heating Distributions, Rei d = 8.23x 1 0 4 
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Figure 9b. Grid Resolution and Grid Adaptation 
Effects on Computed Wake Heating Distributions, 
Re liD = 8.23x104 
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Figure 10. Sample Computed Wake Streamlines, Reip = 8.23xl0 4 
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Figure 11. Reynolds Number Effects on Computed 
Heating Distributions 



Figure 12a. Reynolds Number Effects on Computed 
Heating Distributions - Forebody Details 
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Figure 12b. Reynolds Number Effects on Computed 
Heating Distributions - Wake Details 



Figure 13a. Experiment vs. Computation, 
Forebody Heating, Rei p = 8.23xl0 4 



Figure 13b. Experiment vs. Computation, 
Wake Heating, Re i d = 8.23xl0 4 



Figure 14a. Experiment vs. Computation, 
Forebody Heating, Rej d = 1.52xl0 5 


13 

American Institute of Aeronautics and Astronautics 









Figure 14b. Experiment vs. Computation, 
Wake Heating, Re id = 1.52xl0 5 
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Figure 15a. Experiment vs. Computation, 
Forebody Heating, Rej d = 3.15xl0 5 



Figure 15b. Experiment vs. Computation, 
Wake Heating, Rej d = 3.15xl0 5 
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Figure 17. Computed Streamlines vs. Oil Flow Patterns, 
Rej D = 8.23x104 



x (m) 

Figure 18. Local Knudsen Number Contours, 
Re liD = 8.23x104 



Figure 19. Local Knudsen Number Contours, 
Rej D = 1.62x105 



Figure 20. Local Knudsen Number Contours, 
Re 1;D = 3.15x105 
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